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We propose a powerful method of controlling interaction between silicene and a substrate utilizing 
quantum size effect, which allows to grow silicene with tailored electronic properties. As an example 
we consider silicene on ultrathin Pb(lll) layers, and demonstrate how the properties of silicene, 
including the binding energy, and the Dirac bands, can easily be tuned by quantum well states 
of the substrate. We also discover a novel mechanism of protecting the Dirac electrons from the 
influence of the substrate. This is associated with special arrangement of a part of Si atoms in 
silicene. These findings emphasize the essential role of interfacial coupling and open new routes to 
create silicene-like two-dimensional structures with controlled electronic properties. 
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Silicene, the two-dimensional (2D) allotrope of silicon, 
has attracted considerable attention due to its similar to 
graphene exceptional electronic properties [Ei- How¬ 
ever, silicene and graphene differ in the atomic structure 
- in both of them atoms are ordered in the honeycomb 
lattice, but due to sp^ hybridization, graphene is entirely 
flat, while mixing of sp^ and sp^ hybridization in silicene 
leads to low-buckled structure. Due to this buckled struc¬ 
ture a band gap can be easily generated and controlled 
by electric field H. Properties of silicene can also be 
tuned by chemical functionalization [bl-Q. 


Since silicon tends to hybridize in sp^ form there is no 
graphite-like layered material composed of silicon atoms, 
thereby epitaxy is the only way to obtain the silicene 
layer, as it was successfully demonstrated on a few dif¬ 
ferent substrates lilil. However, presence of the sub¬ 
strate may and usually does alter the properties of the 
free-standing silicene. Indeed, despite extensive studies 


1 , despit 

of silicene on Ag(lll) surface [qL Il4l - l2l| problem of band 
structure (presence or absence of linear bands) does not 
yet have a precise answer. There are firm arguments that 
binding of silicene to Ag surface is so strong that it de¬ 
stroys not only the Dirac cone but also 2D character of 
silicene. However, a recent paper of Y. Feng et al. [22| 
indicates that Dirac cones survive and appear in K points 
of silicene 3x3 Brillouin zone. These studies clearly point 
to the importance of the silicene-substrate interaction. 

Thus searching of alternative substrates for the silicene 
growth is a still debated issue. Semi-conducting surfaces 
seem to be particularly attractive in view of their weak 
interaction with silicene layer, which is beneficial for pre¬ 
serving the linear nature of the band structure 19|, 23- 


l27l |. However, such weak interaction will likely result in 
clustering of silicon atoms, and may prohibit formation 
of 2D Si layer at all. In the case of metal substrates, the 
silicene-substrate interaction is too strong to maintain 
the properties of the free-standing silicene, and leads to 
substantial modifications or to destruction of the Dirac- 


fermion spectrum [28|. It seems that the best choice 
would be the substrate featuring moderate interaction 
with silicene. Recently, we proposed to use lead as a sub¬ 
strate, a material which seems to fulfill the above inter¬ 
action requirement [1^. Indeed, in this case, the binding 
energy falls between energies characteristic of silicene on 
a typical semiconductor and on a metal substrate. As a 
result, the Dirac tt bands are only slightly modified with 
the main contribution of the silicene 3pz orbitals. 

Searching for new templates is one of the most natu¬ 
ral directions in the field of silicene formation. However, 
it is rather laborious, often complicated, and not much 
effective. In the present work, we propose an alternative 
approach, utilizing quantum phenomena rather than new 
materials, to form silicene, and to control its properties. 
The idea is to use as a substrate metallic quantum wells 
(QW), i.e. ultrathin metal layers, in which a quantum 
size effect (QSE) takes place. It is well known that the 
QSE influences the properties of thin metal films and 
makes the quantities related to atomic structure, ener¬ 
getics and electronic properties oscillatory functions of a 
number of monolayers |3Ql435j . 

The positions of QW states are determined by thick¬ 
ness of the slab, which is very attractive in a view of 
silicene growth. It should be possible to manipulate the 
strength of the silicene-QW interaction, thus to allow to 
grow silicene with tailored electronic properties. As an 
example we consider silicene on ultrathin Pb(lll) layers, 
and demonstrate how the properties of silicene, including 
binding energy, and the Dirac cone, can easily be mod¬ 
ified by QW states. The proposed idea can be applied 
to other 2D materials, like germanene, stanene, etc., on 
different QSE substrates. Eurthermore, we also discover 
novel mechanism of protecting the Dirac electrons from 
influence of the substrate, which is associated with the 
presence of Si atoms sticking out of silicene layer. These 
findings emphasize the essential role of interfacial cou¬ 
pling and open new routes to create 2D structures with 
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controlled electronic properties. 

The calculations were performed within density func¬ 
tional theory using SIESTA code sEfioj. The general¬ 
ized gradient approximation in Perdew-Burke-Ernzerhof 
form 0 was utilized as the exchange-correlation func¬ 
tional. The electron-ion interactions were represented 
by the Troullier-Martins norm-conserving pseudopoten¬ 
tial (ill . The plane wave cutoff for all calculations was set 
to 200 Ry. The Brilloiun zone was sampled by 12x12x1 
k points, according to the Monkhorst-Pack scheme (i^ . 

To avoid additional strain we used superlattices 3x3 
of Pb(lll) surface and a/7 x \fl of silicene with approxi¬ 
mately 3% of lattice mismatch. The slab contains differ¬ 
ent number of Pb monolayers (from 2 to 8) and a silicene 
layer on the top. Only atoms from the lowest layer were 
fixed in their bulk positions and the rest of the system 
was fully relaxed until all the forces acting on atoms were 
lower than 0.01 eV/A . The initial geometry of each QW 
system was the corresponding relaxed Pb slab with the 
bulk in-plane lattice constant and the free-standing sil¬ 
icene layer. 

The STM simulations have been performed within the 
Tersoff-Hamann approach in the constant-current mode 
3 EH • The analysis of charge transfer between different 


atoms has been done according to Bader 0. 

Eigure[T](a) (blue line) shows the binding energy of sil¬ 
icene to Pb(lll) film composed of a different number of 
layers. Clear oscillations of as a function of Pb thick¬ 
ness are visible. This is the result of different positions 
of quantized energy levels of Pb film due to the quantum 
size effect (see Eig. 1 in Ref. 0). Obviously, these 
quantum well states are also responsible for binding en¬ 
ergy El of the top Pb layer to rest of Pb slab [see red line 
in Eig. [11(a)]. Thus both, E^ and Ei, show oscillatory 
behavior as a function of number of Pb layers, but they 
oscillate in antiphase. It means that if the top Pb layer 
is bound strongly (weakly) to rest of the slab, it binds 
silicene weakly (strongly). This is an example of common 
wisdom in chemistry on the bonding of atoms, i.e. if a 
given atom is strongly involved in bonding with one of 
its neighbors, at the same time it weakens its bonds with 
the other neighbor EH- Such behavior of the binding 
energy can be exploited in the growth of silicene layer. 
Changing the thickness of a slab, thus the binding energy, 
we modify the growing conditions for silicene. Thus we 
are equipped with a new tool to control the formation of 
silicene. The same scenario should be realized in other 
2D systems on various thin metal layers, provided van 
der Waals interactions are not much important, and sub¬ 
strates exhibit QSE. Obviously, this is unlikely to work 
for graphene. 

Surprisingly, the binding energy oscillations are not re¬ 
flected in substantial changes of the silicene atomic struc¬ 
ture. The silicene layer is buckled and consists of Si 
atoms in their top and bottom positions [IH- The ar¬ 
rangement of Si atoms is almost the same, independent 
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FIG. 1: (Color online) (a) The binding energy of silicene to 
Pb films (blue line), and of top Pb layer to rest of the slab 
(red line) as a function of different number of Pb layers. Evo¬ 
lution of energy gap at the Dirac point (b), position of the 
Dirac point (c), and charge transfer between silicene and Pb 
substrate (d) with Pb QW thickness. 


of QW thickness. The influence of QSE leads only to dif¬ 
ferent mean separations of silicene and Pb surface, with 
large (small) distance for weak (strong) silicene-substrate 
interaction, and to small variations of the buckling, only 
to 0.06 A. The energy ordering of different silicene recon¬ 
structions, predicted in Ref. [29|, is not sensitive to QSE 
either. The lack of structural modifications is supported 
by recent experiment on silicene on Ag(lll) films, albeit 
containing more atomic layers EH- 

Scanning tunneling microscopy (STM) topography im¬ 
ages calculated for silicene on quantum wells of different 
thickness show, however, noticeable changes, as one can 
deduce from experiments on bare Pb layers EH- These 
changes are most pronounced for Si atoms in the top 
position, as it is visible in Eig. [2] for marked bright pro¬ 
trusion. This suggests that electronic effects play an im¬ 
portant role in the system. Namely, the QW states of 
the substrate provoke changes in the electronic structure 
of silicene, thus influence the STM topography. Indeed, 
this is clearly visible in Eig. |3l where evolution of the en¬ 
ergy bands as a function of Pb thickness is shown. In all 
cases bands with linear dispersion can be detected, albeit 
modified. Due to the presence of a substrate, the inver¬ 
sion symmetry of silicene lattice is broken, which results 
in the opening of an energy gap at the Dirac point. The 
magnitude of this gap, is related to QSE via the binding 
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FIG. 2: (Color online) (a) Constant current STM topography of silicene on 3 ML thick Pb film, calculated at the bias voltage 
-IV. (b) and (c) corresponding images for 4 and 5 ML. Arrows point to the bright protrusions, most strongly affected by QSE. 



FIG. 3: Evolution of band structure silicene on Pb layers as a function of Pb thickness. Black solid lines show the unsupported 
silicene band structure, which has been calculated by fixing Si atoms positions as obtained in full slab calculations and removing 
the substrate. 


energy and the mean silicene-Pb surface separation, and 
oscillates as a function of QW thickness [see Fig. [I](b)]. 

Similarly, position of the Dirac point [Fig. [11(c)], de¬ 
fined as energy in the middle of the Dirac cone gap, 
can be shifted by approximately 0.4 eV by varying QW 
thickness. This effect is accompanied by charge trans¬ 
fer between both subsystems, as it is evidenced in Fig. 
H] (d). However, the situation is more complicated, and 
additional mechanism must be present, since this charge 
transfer is too low to shift the Dirac point so deep in 
the valence bands. Indeed, not all Si atoms receive neg¬ 
ative charge, as the Bader analysis indicates. Namely, 
Si atoms in the lower layer, which contribute mainly to 
the Dirac cone, are doped by electrons from the top Si 
layer, which further shifts the Dirac cone down. Thus the 


electron charge is transferred from both Pb and top Si 
atoms to the bottom Si atoms. Moreover, there is strong 
hybridization of electron states, mainly derived from p- 
orbitals of Si (a and tt bands of silicene) and Pb atoms, in 
vicinity of the F point. This can be observed while com¬ 
paring unsupported silicene bands with the band struc¬ 
ture of the full slab near the F point. Clearly, strong 
modifications of the unsupported silicene bands are visi¬ 
ble. Note that the strain imposed on silicene layer is too 
weak to shift the Dirac cone [5l| . 

The charge transfer between top and bottom Si atoms 
not only shifts the Dirac cone, but is an intrinsic fac¬ 
tor of a novel effect, which we call self-protection of 
silicene. The important role in the mechanism of self¬ 
protection falls to the top Si atoms in supported silicene 
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layer. Namely, the tt electrons of silicene try to avoid the 
destructive influence of a substrate as much as they can 
for the price of strong interaction between top Si atoms 
with Pb film. Such behavior is triggered by the QSE, 
as the electron wave function forms a standing wave lim¬ 
ited by bottom of the slab and the top-most objects (top 
Si atoms). Thus the role of the top Si atoms is to de¬ 
couple the Dirac electrons in silicene from other parts 
of the system. Consequently, we can think of silicene 
on a substrate as an electron doped quasi-free-standing 
layer with top Si atoms acting as substitutional impuri¬ 
ties. This scenario is supported by projected band struc¬ 
ture, shown in Fig. IH Clearly, the main contribution 
to the linear bands comes from the bottom Si atoms, 
with a weak impact of Pb 6pz and even weaker of top 
Si 3pz states. On the other hand, there is strong inter¬ 
action between Pb and top Si Pz orbitals, forming QW 
states. This is a hallmark of the self-protection mecha¬ 
nism. Note that there is no trace of this QW state in the 
electronic structure of the bottom Si atoms. Thus the 
Dirac cone in silicene on QW survives and the tt elec¬ 
trons of silicene try to be decoupled from the rest of the 
system owing to strong interaction between Pb and top 
Si atoms. We think that the self-protection scenario is 
realized in other supported 2D systems, like silicene on 
Ag(lll), where the atomic structure of silicene also fea¬ 
tures a few Si atoms in the unit cell sticking out of the 
Si layer. However in this case the self-protection mecha¬ 
nism may be too weak to compete with the very strong 
silicene-substrate interaction. 

The decoupling of 7r-electrons depends on relative po¬ 
sition of the Dirac point and QW state, and can be 
slightly suppressed if both coincide, as it is seen in Fig. 
[5l This may occasionally happen for certain thicknesses, 
and rather for thicker films, where the separation between 
QW states is smaller. 

Finally, we would like to discuss the experimental re¬ 
alization of the proposed method of growing and con¬ 
trolling 2D silicene-like systems. The strongest effect 
should be observed for thin films. However, the QWs 
also require substrates. Therefore, it is important to 
chose the substrate on which the film can be grown in 
a layer by layer mode. Furthermore, present calculations 
have been performed for unsupported QWs. The un¬ 
derlying substrate on which thin films could be grown 
will certainly influence the QW states to some extent. 
Therefore, one can choose a substrate with minor impact 
on QW states like, for example, single-crystalline epi¬ 
taxial graphite Q. On the other hand, the supported 
quantum wells may even help to control the properties 
of silicene, as in the case of Pb films on Si(lll) surface, 
where the QW states become less dispersive, and remain 
flat in a large part of the Brillouin zone [35[. Thus the 
silicene properties should be more controllable. Finally, 
we would like to stress that the proposed idea can be 
utilized in other 2D materials on various substrates ex¬ 


hibiting QSE. 

In conclusion, we proposed a powerful method of con¬ 
trolling the silicene-substrate interaction owing to the 
quantum size effect. We have demonstrated that the elec¬ 
tronic properties of silicene on metallic quantum wells, 
like binding energy, position of the Dirac cone and mag¬ 
nitude of the energy gap at the Dirac point, can easily be 
tuned by quantum well states of the substrate. This idea 
can be utilized in other 2D silicene-like materials grow¬ 
ing on various thin films. We also discovered the effect 
of self-protection of silicene, in which the top Si atoms 
play a crucial role, taking on the interaction with the 
substrate, and decoupling Dirac electrons from the rest 
of the system. These findings emphasize the essential 
role of interfacial coupling and open new routes to create 
two-dimesnional structures with tailored electronic prop¬ 
erties. 
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FIG. 4: (Color online) Projected band structure of silicene on 4 ML thick Pb film, (a) Contribution of Pb 6pz, (b) top-Si Sp^, 
and (c) bottom-Si Sp^ states to the electronic bands. QWS and DC stand for the quantum well state and the Dirac cone, 
respectively. 



FIG. 5: (Color online) The same as in Fig. |4]for silicene on 5 ML thick Pb film. 
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